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1  .  INTRODUCTION  j 

^  The  positioning  of  nulls  in  an  antenna  array  field  pattern  is  essential  to 
the  performance  of  the  antenna,  in  being  capable  of  blocking  interference. 
The  null  placement  must  be  acnieved  in  such  a  way  that  the  field  pattern 
in  other  directions  is  not  adversely  affected. 

One  of  the  most  efficient  methods  of  null  placement  is  by  perturbing  only 
the  phases  of  the  array  elements.  iiaP6)— presentiytwo  approaches  to  the 
placement  of  nulls  by  phase  perturbation.  The  first  is  a  least  squares 
method  based  on  exact  or  approximate  null  placement,  applicable  to  one¬ 
dimensional  arrays  and  extendable  to  two-dimensional  arrays,  developed  for 
real  quiescent  patterns  which  apparently  allows  polygonal  arrays  (in  this 
study,  octagonal  arrays)  to  be  considered.  The  second  is  a  minimax  method 
in  one  or  two  dimensions  based  on  null  placement,  which  readily  permits 
the  omission  of  failed  elements  and  which  involves  only  the  perturbation 
of  selected  element  phases  or  amplitudes.  ^  u  ■  j>  '  '■  '  ^  , 

We  believe  that  these  two  app ’caches  can  provide  between  them  a  versatile 
choice.  The  least  squares  method  is  extremely  efficient,  the  number  of 
parameters  being  simply  the  number  of  null  constraints,  but  ail  antenna 
weights  need  to  be  perturbed.  The  minimax  method  is  much  more  expensive, 
since  all  antenna  weights  are  parameters,  but  often  only  a  small  number  of 
weights  need  to  be  changed.  Both  approaches  are  fairly  robust.  The  least 
squares  approach  apparently  permits  polygonal  arrays  to  be  adopted,  and 
probably  more  general  geometries  and  configurations  of  failed  elements 
might  be  adopted.  The  minimax  approach  is  based  on  a  very  general 
optimisation  algorithm  and  therefore  in  principle  permits  rather  wide 
ranging  specifications  of  constraints  to  be  imposed  by  the  user. 

Both  approaches  are  still  under  develooment;  the  theory  is  incomplete  and 
algorithms  have  not  been  fully  tested.  In  particular  r.o  a  priori 
guarantee  can  be  given  that  the  least  squares  technique  will  work  for  any 
given  polygonal  array,  and  at  present  it  cannot  be  guaranteed  that  the 
minimax  algorithm  will  always  converge  or  that  it  will  lead  to  a  minimal 
number  of  phase  changes.  Since  the  current  AFOSR  research  grant  is  not  to 
continue  for  a  second  year,  as  originally  planned,  it  will  not  be  possible 
to  complete  this  work  under  AFOSR  support  during  the  coming  year. 


2.  PHASE-ONLY  NULLING  BY  LEAST  SQUARES  METHODS  IN  ONE  DIMENSION 

We  consider  a  linear  array  of  N  isotropic  elements  as  shown  in  figure 
The  field  pattern  is  given  by, 


N  id  u 
P0(u)  =  I  V  *  . 
n  =  1 

tn 

The  weight  a^  is  the  complex  excitation  of  the  n  element.  The  phase 
reference  is  taken  to  be  the  centre  of  the  array,  hence  the  weights  d 
are  given  by, 


and 


d  =  (N-l)  -(n-1 )  -  -d  . ,  n=1 ,  .  . N 

n  — N-n+1 

u  =  2ird  sine, 

X 


where  X 

d 

9 

The  interelement  spacing 


«  wavelength, 

=  interelement  spacing, 

=  angle  subtended  with  the  normal  to  the  arra 
is  taken  to  be  half  the  wavelength  throughout 

d 

r  <  »| 

•  i 
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Figure  1  Geometry  of  the  Array. 


We  require  a  set,  $  ,  n=i,..N  of  phase  perturbations  which  impose  nulls  in 
required  locations  whilst  retaining  the  characteristics  of  the  quiescent 
pattern.  Let  u  ,  k=1,..K,  be  the  directions  in  which  nulls  are  required, 
then  to  obtain  $  ,  n=1..M,  we  minimize  the  integral  of  the  square  of  the 
discrepancy  between  the  perturbed  pattern  and  the  original  pattern,  which 
is  readily  shown  to  equal 

N  i* 

F  =  l  c  a  (e  n  -  I)2,  (2.1) 

n*1 

(c  are  assumed  to  be  real  and  positive).  To  impose  the  nulls  in  the 
required  locations  the  phase  perturbations  must  satisfy  the  following 
constraints , 


i$ 


n= 


a  e 
,  n 


n  idnuk 
e 


0,  k=l,..K. 


(2.2) 


The  problem  is  nonlinear  in  general  and  cannot  be  solved  analytically,  but 
numerical  solutions  can  be  obtained  using  nonlinear  programming. 

We  consider  the  case  when  the  quiescent  pattern  is  real,  and  hence  the 
element  weights  are  conjugate  symmetric, 

aN-n*1  =  an'  n=1 

in  this  case  it  can  be  shown  (Shore  [1984a])  that  the  required  phase 
perturbations  are  odd-symmetri^ . 


N-n+1 


n 


(2.3) 


Writing  the  coefficients  in  the  form 


a 

n 


'idnus 


then  due  to  the  odd  symmetry  of  the  pertubations  $  ,  and  the  coefficients 
d  ,  n=l,..N,  the  constraints  given  in  equation  (2.2)  can  be  written. 


N' 

C.  -  ~  la  I  cost*  +  d  (u.  -  u  ) } ,  k*1 , . .X. 
k  “ .  1  n 1  n  n  k  s 

n= 

(2.4) 


The  objective  function  can  be  rewritten  in  the  form, 


N 

F  =  7  c  la  1 2 ( 1  -cos<s  ) 
L ,  n  1  n 1  n 

n=  1 


(2.5) 


Given  the  objective  function  ?,  equation  (2.5)  and  the  cc 
equation  (2.4),  we  define  the  Lagrangian  L  as  follows, 


■.straints  , 


K 


L  =  F 


(2.6) 


where  the  coefficients  A ,  k=1,..K  are  (real)  Lagrangian  multipliers.  A 
necessary  condition  for  the  perturbations  <p  to  locally  minimize  F  subject 
to  the  constraints  C^=0,  is  that  the  partial  derivatives  of  L  with  respect 
to  the  (p^  be  zero.  Hence  a  necessary  condition  for  there  to  be  a  minimum 
is  that 


K 

9F  -  y  A  3C,  =  0 ,  n  =  1  , .  .N. 
3*  k-1  3^ 

n  1  n 


(2.7) 


Differentiating  F 


and  with  respect  to  $  , 


say,  gives, 


2c 


pi  “pi 


"  "rP* 


n  — I 

f  1 


%» 


(2.8) 


and 


(2.9) 


lap  |sin[<pp  *  dp(u^  -  us)] 

p  =  1 , . .N, 
k=*1  ,  .  .K. 


Substituting  equations  (2.8)  and  (2.9)  int">  (2.7)  gives. 


K 

2ep  I  ap  I  2  sin*p  =  '  l  \k  I  ap  I  sidL'Pp  +  a?(Jk'Js'  - 
k«1 


A  little  algebraic  manipulation  gives  <+>_  in  the  form 

K 


K 

tan  <fp  =  -  E^  Aksin[dp( jk-u2 )  ]  (2.10) 

K 

2cn I an !  +  :  xk'Os[dp(uk-u£)] 

k-1 

Comoaring  this  form  for  the  phase  perturbations  <j>  with  that  given  by 

P 

Shore  [1983],  it  is  clear  that  the  coefficients  which  Shore  refers  to  as 
the  'beam  coefficients'  are,  in  this  case,  the  negative  Lagrangian 
multipliers.  (In  Shore  because  of  the  slightly  different  form  given  for 
the  $  ,  n=1,..N,  the  coefficient;  are  actually  the  negative  Lagrangian 
multipliers  divided  by  two.)  This  is  the  case  only  when  the  pattern  is 
real,  and  hence  we  can  take  advantage  of  the  symmetry  in  the  coefficients 
ensuring  the  constraints  are  real. 

The  unknowns  in  the  minimization  problem  are  X,  ,  k=1,..K  ,  the  problem 
size  is  dependent  upon  the  number  of  constraints  K,  rather  than  the  number 
of  array  elements  N.  This  is  clearly  advantageous  as  the  nui?K°"  of 
elements  is  generally  much  larger  than  the  number  of  constraints  required. 


At  present  the  problem  is  solves  by  using  a  WAG  routine  (ZO^VDF)  wnicr. 
employs  a  sequential  quadratic  programming  algorithm  (3QP)  (see  Gill, 
Murray  and  Wright  [1981]).  The  routine  requires  an  initial  estimate  for 
the  solution  and  routines  to  evaluate  the  objective  function,  constraints 
and  their  derivatives  with  respect  to  the  coef f icie.nts . 

The  results  given  in  table  1  are  for  a  problem  in  which  there  are  -■ 

antenna  elements,  c  =  la  I  =  1,  n=1..N,  and  the  quiescent  beam  direction 
n  1  n‘ 

u  =0.  The  null  is  located  at  0,  =0 . ^3^33  (rads)  and  the  tolere.nce 

s  K 

levels  required  by  the  routine  were  set  to  le-10.  This  example  is  taxer, 
from  Shore  [1983],  and  as  stated  previously  the  'beam  coefficients'  in 
the  formulation  presented  here  are  twice  those  of  Shore  and  opposite  :n 
sign,  but  the  resulting  weight  perturbations  are  the  same. 


2.1  Linearization 

If  the  phase  perturbations  are  assumed  small,  then  we  can  employ  the 
following  approximations, 


tan (<f)  ■=•  <p, 


and 


i  <Pr 


1  +  if 


The  weight  perturbations  w  =*  a  (e '  ^  can  then  be  approximated 


n  n 


wn  =  ian*n 


(2.11) 


■* 


This  form  for  the  weights  allows  us  to  rewrite  toe  cancellation  pattern, 
in  terms  of  the  new  coefficients,  as  the  sum  of  two  beams 


Ap(u)  =  \  l  X,  V  J_Teldn(u  ^2us  Jk  J ;  -  eidn(u  uk  5 
k  =  1  n  =  1  cn'- 


(2.12) 


where , 


c’  =  2cn  -  1  >  ^cos[dn(uk-us)] 

I an |  k=i 


One  beam  is  in  the  direction  of  the  required  null,  which  cancels  the 
original  pattern,  and  the  second  is  in  the  symmetrical  location  with 
respect  to  the  main  beam;  which  leads  to  an  enhancement  of  the  pattern  at 
this  point.  Pattern  enhancement  occurs  in  the  example  described  above  an 
is  illustrated  in  figure  2,  tr.e  null  at  uk  =  0.^3633  (rads)  is  reflected 
by  pattern  enhancement  at  u  =  -0.^3633  (rads). 


It  should  be  noted  that  if  we  can  make  the  approximations  for  #  and  the 
weights  as  shown  above,  and  if  the  coefficients  Ak  are  small  relative  to 
2,  with  cn  =  |an|  =  1,  then  by  neglecting  any  contributions  from  other 
cancellation  beams  at  u=uk ,  we  can  obtain  an  estimate  for  the  coefficient 
as  follows;  the  cancellation  pattern  can  be  approximated  by 

Ap(u,  )  »  -  -  A,  ,  k  =  l  ,  .  .K 
K  14  K 

and  since  by  definition  wo  have 
Ap(uk)  =  -  pQ(uk) 

then  the  coefficients  can  be  approximated  by 

Ak  =  Up0(uk) . 

N 


(2.13) 


Table  2a  shows  the  results  wher.  nulls  are  placed  as  6  =  0.59556,  0.65575, 
0.71887  (rad),  the  approximate  coefficients  in  this  case  are  -0.12646, 
0.11923,  ana  -0.11349  respectively,  which  are  clearly  of  the  same  order  as 
the  calculated  values.  However,  table  2b  illustrates  the  effect  of 
imposing  beams  close  together.  As  locations  of  nulls  appro  ch  each  other 
they  have  a  greater  influence  on  each  other  and  the  coefficients  are  no 
longer  objectively  independent.  This  is  also  true  when  imposing 
additional  nulls  close  to  locations  symmetrical  to  the  original  nulls. 


2.2.  Symmetrical  Nulls 


It  is  not  possible  to  synthesize  nulls  at  symmetrical  locations  using  the 
linearized  form  and  approximation  to  the  weights.  (A  proof  is  included  in 
appendix  1  for  completeness).  When  placing  symmetrical  nulls  no 
assumption  can  be  made  about  the  size  of  the  phase  perturbations.  Shore 
[1984b]  considers  the  problem  of  symmetrical  nulls,  and  concludes  that  it 
is  always  possible  to  achieve  symmetrical  nulls  in  a  linear  array  pattern, 
but  it  may  result  in  interference  patterns.  We  expand  these  ideas  here 
and  show  that  in  some  cases  it  is  not  possible  to  place  symmetrical  nulls 
using  this  method. 

If  we  consider  the  constraints  and  the  for  symmetrical  null  locations 
uk  =  ±u»,  assuming  ug  =  0,  we  have  from  equations  (2.4)  and  (2.10); 

N 

C,  -  l  cos(4>n  +  dnu*)  =  0, 
n=1 


(2.14) 


and 


N 

C 2  "  I  cos ($n  -  dnu*)  =  0, 
n=1 


tan4>n  = 


sin(dnu*) ( A1 -A?) 


2cn I an 1 


+  cos(dru*)  (  X1  a ^ ) 


(2.15) 


Using  equations  (2.14)  and  (2.15)  we  would  like  to  find  the  conditions 
which  specify  whether  there  are  an  infinite  number  of  solutions,  no 
solutions  or  a  unique  solution.  Shore  [1984b]  states  that,  provided  that 
the  phase  perturbations  are  not  restricted  to  be  small,  nulls  can  be 
imposed  at  locations  symmetrical  about  the  main  team;  this  is  not  always 
the  case  and  can  easily  be  proven.  Taking  the  case  when  cn  =  |an|  =  1 
with  0  =  ir/6  (rads),  u*  =  irsinQ  »  ir/2.  From  equation  (2.15)  we  have 


When  dn 
when  dn 


tan 


sin(dntt/2)  ( 1 ,  -X 2 ) _ 

2  +  cos(dn/2)(A1+X2) 


n=1 , . .N. 


is  even,  sin  (dnTt/2j  =  0,  and  therefore  $n 
is  odd,  cos(dnir/2)  *  0,  and  therefore  <f>n  = 


(2.16) 


t tt ,  t  is  an  integer; 
±  tan-1  UrX2}/2. 


The  terms  in  the  summation  of  the  constraints  given  by  equation  (2.15) 
are : 

(i)  for  even  dn,  cos  ( <|!n±dnif/2)  =  cosCtjir)  =  ±  1, 

(ii)  for  odd  dn,  cos(<fin±dn7r/2)  =  cos(T±tt/2), 
where  T  =  tan  ^(Ai_A2)/2.  The  summation  of  the  terms  results  in  the  pair 
of  equations 


^li.  COsCTr— )  *  '<  =  0  , 

2  2 

where  k  is  a  non-zero  integer.  This  leads  to  the  equations, 

s^nT  =  ±  -  , 

N-i 

which  indicates  that  the  cor.strair.cs  are  two  parallel  lines  and  thus  can 
never  both  be  satisfied.  In  this  case  there  is  no  solution  to  the 
constrained  minimization  problem.  (Here  it  is  assumed  that  (N-1)/2  is 
odd,  however  the  result  is  of  the  same  form  if  (N-l)/2  is  even). 

In  general  there  is  a  solution  to  the  constrained  minimization  problem  and 
the  amount  of  interference  which  occurs  is  dependent  on  the  position  of 
the  null  relative  to  a  null  in  the  quiescent  pattern.  Figure  3 
illustrates  this  relationship  by  placing  nulls  at  intervals  between 
quiescent  nulls  of  a  41  element  array. 

There  appears  to  be  a  relationship  between  the  optimal  coefficients  and 
the  unperturbed  pattern  value  but  we  have  been  unable  to  establish  that 
relationship  at  present. 


2.3-  Increasing  null  width  by  Higher  Order  Constraints 

Since  the  placing  of  nulls  at  locations  close  together  results  in 
cancellation  beams  interference,  as  illustrated  in  the  previous  section, 
it  may  prove  wise  to  employ  alternative  methods  for  increasing  the  width 
of  a  null. 


pa(uk)  -  0, 


k=1 , .  .K 
u=0 ,  .  .M 


(2.17) 


wnere  uk  denotes  tne  _ocatior.  tr.e  .\  intent  erer.C'3  a:rec:;:r.s. 

Previously  we  have  considered  v=0,  however,  by  including  higner  crier 
derivatives,  the  null  width  is  broadened.  This  was  illustrated  for  null 
synthesis  with  phase  end  amplitude  perturbations  by  Steyskal  [1982]. 

Here,  we  investigate  the  use  of  higher  order  constraints  in  the  context  of 
phase-only  nulling. 


Clearly  now  the  total  number  of  constraints  is  P-K(>'.+  1 )  ,  and  they  are 
given  by 


C  =  l  (dn)v  |  an  i  cos[4r_  *  dn  (uk-us)j 


(2.18) 


for  even  v, 


cp  =  l  (dn)V  laJ  sir-  C*n  +  dn  (uk“us)] 


(2.19) 


for  odd  u. 


The  Lagrangian  is  now  of  the  form 
P 

L  ■  F  ~  l  *pCp  , 

p  =  1 


and  again  the  condition  for  d  minimum  is 

-  l  xp  l£p  -  0 
P”'1  34>n 


(2.20) 


From  equations  (2.18)  and  (2.19)  the  derivative  of  the  constraints  are 
given  by:- 


-  |aj  (dnr  sin  Un  ♦  dn(uk-us)] 


(2.21) 


for  even  v, 


|an|(dn)V  cos  Un  -  dn(uk~Js)] 


(2.22) 


for  odd  v. 


From  equations  (2.9),  (2.20),  (2.21)  and  (2.22)  we  have 

K  M 

2cnlanlsin*n  “  “E  l  xkv(dn)V  sin  [^n  +  dn(uk"us)] 

k=1  v=0 

step  2 


K  M 

+  £  l  Xkv;(dn)V  cos  C*n+dn(uk"us)]’ 

k=1  v=1 

step  2 
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where  kv  =»  (k-1 )  (M+l  )  *•(  u*!) .  After  a  little  algebraic  manipulation  we 
obtain, 

K  M 

^  ^k\Zdn^  3indn^  Jn"Js  ^  4  4  xku^dn^  cosdn(  Jk'u5 ) 

k=l  v=0,2..  k  =  1  v=*1  ,  3.  . 

step  2  step  2 


tan  <j> 


n 


KM  K  M 

2cjan|  +  l  l  Akv(dn)vcosdri(ui<-u2)  -  I  I  V/ ^ ) 'Jsindn  (uk-us  ) 


k  =  1  v=0 ,  2 . 


k=1  v= ' , 3 • 


(2.23) 


If  v  =  0  then  (2.23)  reduces  to  equation  (2.7). 

The  problem  can  be  solved  in  the  same  manner  as  previously  except  now 
there  are  K(M+1)  variables  instead  of  K.  The  method  used  to  solve  the 
constrained  minimization  problem  requires  derivatives  of  both  the 
objective  function  and  the  constraints,  the  equations  for  these  are  given 
in  appendix  II. 


2.3.1  Numerical  Results 

The  effect  of  including  higher  order  constraints  is  illustrated  in  figures 

*»,  the  corresponding  coefficients  are  given  In  table  3.  Fig-are  4(a) 

illustrates  the  quiescent  pattern,  with  31  elements  and  u  =0  and 

s 

an»1  n»1,...N.  Figure  4(b)  illustrates  the  pattern  with  a  zero  order  null 
located  at  uk«0.3,  the  null  is  indicated  by  the  vertical  line.  Figure 
4(c)  and  (d)  illustrate  clearly  how  the  null  is  broadened  by  the  addition 
of  higher  order  nulls,  with  a  first  order  and  second  order  null 
illustrated  respectively.  It  is  evident  from  the  coefficients  in  Table  3 
tha«-  in  this  example  the  original  beam  coefficient  does  not  vary  greatly 
with  the  addition  of  higher  order  constraints.  However  the  whole  pattern 
is  affected  by  the  higher  order  constraint,  resulting  in  some  amount  of 
pattern  enhancement  on  the  half  range  symmetrical  to  the  null  location. 


3. 


PHASE  ONLY  NULLING  BY  LEAST  SQUARES  METHODS  IN  TWO  DIMENSIONS 


3 . 1  Notation 


For  a  two-dimensional  planar  array  the  phase  equation  is  of  the  form, 


ek,£  *  —  dx  wxk  u  +  —  dy  wy2.  v 
X  x 


(3.1  ) 


k= 1 . .Ny,  JU1 , . .Nv, 

x  y 


and  the  resulting  field  pattern  is  given  by 


P  = 


Y 

L 


n=1  m=1 


an  ,m 


gi9n,m 


(3.2) 


Here 


N 

N 


x 

y 


wxk 


w  Vi 

k,l  or 


dx’dy 


n,m 


u,v 


=  number  of  elements  in  the  x-direction 
=  number  of  elements  in  the  y-direction 
=  k-(nx+1)/2,  k-1,..Nx 
-  I- (Ny+ 1  )/2 ,  1-1,. .N 

m,n  =  the  reference  coordinates  in  the  x,y  plane 

=  element  spacing  in  the  x,y  directions  respectively 
=  complex  excitation  of  the  n,mth  element 
=  positions  in  the  coordinate  system  shown  below. 


and  u  »  sin  ijj  cos  <p ,  v  =  sin  sin 


where  the  angles  iJj  and  $  are  as  illustrated  in  figure  5. 


We  assume  that  the  element  spaeings  d  and  d  are  half  wavelength. 

x  y 

We  also  assume  that  the  quiescent  pattern  is  real,  and  in  that  case 
the  coefficients  an  m  exhibit  a  symmetry  about  the  centre  element  of 
the  array.  We  shall  return  to  this  point  later. 


3.2  A  Two-Dimensional  Planar  Array 

We  require  a  phase  change  in  the  element  coefficients  which  will 
result  in  the  placement  of  a  null  in  a  given  direction  whilst 
replicating  the  remaining  field  pattern. 

Given  a  set  of  points  {uk,Vk}  k=1,..K  at  which  nulls  are  to  be 
placed,  then  we  require 

p(Uk,  Vk)  =0  k-1 ,..K  (3.3) 

Let  us  consider  the  geometry  of  the  array.  If,  given  a  grid  of 
elements,  the  main  beam  is  assumed  to  be  central  to  the  array,  (that 
is  the  weights  wxn  and  wym  are  as  defined  above),  then  the  phase 
equation  at  each  element  is  as  illustrated  in  figure  6.  (The  example 
here  is  for  a  5x5  grid,  but  clearly  the  basic  pattern  is  the  same  for 
a  general  grid). 

If  the  elements  are  numbered  in  the  order  illustrated  in  figure  6,  it 
can  easily  be  shown  that  in  order  for  the  field  pattern  to  be  real, 
the  complex  coefficients  am  n  must  be  conjugate  symmetric  about  the 
central  element  (in  this  case  element  number  13). 

Ordering  the  coefficients  in  this  manner,  as  a  vector,  we  have 

am,n  *  J  =  m(n-1)  +  m  (3. 4) 

N  =  Nx  x  Ny 


We  can  now  consider  the  field  pattern,  given  in  (3-2)  in  a  vector 

form , 


v  ^  ®  j 

p  =  L  e  J 


(3.5) 


and  in  effect  we  now  have  a  one-dimensional  problem. 


If  the  unknown  phase  perturbat ions  required  to  place  a  null  (or  set 
of  nulls)  is  denoted  by  {$j}  j=1,..N,  then  from  equation  (3.3),  we 
have  the  null  constraints. 


r  i  4>  i  9  k 

)  aj  e  J  e 

j  =  1 


k=1  , .  .K 


(3.6) 


To  ensure  that  the  perturbed  pattern  replicates  the  quiescent  pattern 
everywhere  but  at  the  null  location,  we  minimize  the  sum  of  squares 
of  the  absolute  element  perturbations 


F  -  I  Cj  |  a  j  ( e  ~ 
J  =  1 


(3.7) 


where  Cj  j=1,..N  are  positive  weights. 


Owing  to  the  conjugate  symmetry  of  the  coefficients,  this  can  be 


written  in  the  form, 


F  =  2  l  Cj  |ai  1 2  (1-cos9j  ) , 
j  =  1 


(3.8  j 


The  null  constraints  given  in  (3-6)  can  also  be  re-written  as 


M 

r*  i  o  \.  i  3  ^ 

0  =  Ck  =  l  3j  i  k  2.  j 

j=1 

n 

=  l  |  a  j  j  cos  [  $j  +  9,k]  k=  1  , . .  K  (3.9) 

j-1 

As  in  the  one-dimensional  problem,  we  form  the  La6raiigiau 

K 

L  =  F  -  l  Xp  Cp,  (3.10) 

P-1 

and  the  problem  has  a  'beam  space'  solution  just  as  for  the  one- 
dimensional  problem.  F^om  (3.8),  (3.9)  and  (3-10)  we  can  obtain  the 
relationship 

K 

-  I  Ak  sin  [d  (Uk-Us)] 
k=  1 

tan  (♦  )  =  -  (3.11) 

K 

2  CnKI  +  l  Xu  cos  [d  (Uk-us)] 

k=1 

Details  of  the  algebra  behind  this  relationship  are  as  discussed  in 

§2. 

At  this  stage  it  is  a  simple  matter  to  show  that  all  the  results  for 
the  one-dimensional  case  hold  in  two  dimensions. 

Clearly  the  number  of  variables  in  the  optimization  problem,  namely 

that  of  minimizing  (3.8)  whilst  satisfying  the  constraints  given  in 

(3.9),  is  reduced  from  the  number  of  elements  N  xNv  to  the  number  of 

x  y 

constraints  K,  by  utilizing  (3.11). 
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Just  as  for  the  one-dimensional  array,  in  some  cases  there  is  an 
enhancement  of  the  pattern  in  a  direction  symmetric  to  the  location 
of  the  null.  Examples  of  this  phenomenon  will  be  illustrated  in  the 
numerical  results  given  below. 


3 . 3  Polygonal  arrays 


Billam  (1985)  poses  the  question  of  the  suitability  of  phase  only 
nulling  for  an  octagonal  array  of  elements. 

It  is  possible  to  embed  an  octagonal  array  of  elements  in  a 
rectangular  grid;  as  shown  in  figure  7.  All  the  elements  of  the 
rectangular  array  which  lie  outside  the  octagon  are  clearly  in 
symmetric  positions  about  the  central  clement.  By  putting  the 
initial  weights  of  these  elements  to  zero  and  ensuring  that  they  are 
eliminated  from  further  calculations,  it  is  possible  to  simulate  the 
problem  of  null  placement  in  an  octagonal  antenna  field  pattern. 

This  method  of  embedding  can  in  principle  be  extended  to  any 
polygonal  geometry  of  array. 

We  have  developed  computer  routines  for  the  embedding  of  an  octagonal 
array  into  a  rectangular  array,  and  the  placement  of  nulls  by  phase 
perturbations  to  the  octagonal  array  elements. 

A  listing  of  the  embedding  routine,  which  embeds  the  octagonal  array 
into  a  suitable  rectangle  by  setting  the  appropriate  weights  to  zero, 
can  be  found  in  appendix  III;  it  is  quite  self-explanatory.  In  the 
section  on  numerical  results  below,  we  illustrate  the  difference  in 
the  quiescent  field  for  the  octagonal  and  the  rectangular  arrays, 
show  how  single  and  multiple  nulls  can  be  achieved  for  each,  and  the 
effect  on  other  areas  of  the  pattern. 


3.4  Numerical  Results 


For  the  following  results,  the  optimization  routine  used  to  solve  the 
constrained  minimization  problem  was  taken  from  the  Harwell  library 
of  optimization  routines  ( VF 1 3 ) .  The  method  is  based  on  a  quadratic 
programming  technique  and  is  described  in  Powell  (1982)  and 
Chamberlain  et  al.  (1982).  Linear  approximations  are  made  to  the 
non-linear  constraints,  and  hence  the  placement  of  nulls  at  symmetric 
locations  would  not  be  posible  using  the  routine.  (Tec  §2.2  above). 

The  routine  requires  the  evaluation  of  the  objective  function  and 
constraints,  plus  their  first  derivatives. 

The  examples  below  illustrate  single  and  multiple  null  placement  for 
a  rectangular  antenna  array  and  an  octagonal  antenna  array.  Figures 
8(a)  and  8(b)  illustrate  the  quiescent  sine  pattern  for  both  the 
rectangular  array  and  the  octagonal  array  respectively  with  a  grid  of 
13x13  elements.  The  octagonal  array,  which  is  embedded  into  the 
13x13  grid  of  elements,  has  5  elements  along  each  face,  and  clearly 
the  resulting  pattern  is  more  circular  in  shape. 

For  each  case  we  have  placed  a  null  at  u  =  0.28,  v  =  0.32,  with  a 
—  8 

tolerence  of  10  allowed  on  the  constraint,  and  an  initial  estimate 
of  the  beam  space  coefficient  taken  as  0.1.  Figures  9  and  10 
illustrate  the  resulting  perturbed  patterns  and  the  difference 
between  the  perturbed  and  the  quiescent  pattern  for  a  rectangular 
array  and  an  octagonal  array  respectively.  The  coefficients  for  the 
perturbed  pattern  are  given  in  each  case  in  tables  4  and  5, 
respectively,  (only  the  first  half  of  the  coefficients  need  be  given, 
owing  to  symmetry).  It  is  clear  from  these  tables  that  a  null  placed 
in  the  octagonal  pattern  results  in  a  larger  absolute  beam 
coefficient,  and  this  in  turn  results  in  a  higher  average 
perturbation  but  leads  to  a  deeper  null. 

On  inspecting  the  graphical  results,  it  is  clear  from  the  difference 
patterns  (i.e.  the  differences  between  perturbed  and  quiescent 
patterns)  that  the  octagonal  pattern  is  affected  slightly  more 
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The  least-squares  method  applied  to  both  rectangular  and  octagonal 
array  phase-only  nulling  problems  has  always  given  good  results  in 
the  cases  considered  so  far. 


4.  NULL  PLACEMENT  BY  MINIMAX  METHODS 

Given  an  initial  far  field  pattern 


p0(u) 


idra 


n=  1 


(4.1  ) 


the  perturbed  pattern  becomes 

N 

<-  id„u 

p ( u )  -  l  x„  e  n  '4.2) 

n=  1 

Denoting  the  discrepancy  by  e(u)  then 

M 

e(u)  =  R  (u)  ♦  i  I  ( u )  =  l  (an-xn)eidr,U  (4.3) 

n=  1 


For  amplitude-only  perturbation  xn/an  is  real  (and  so  is  an)t  and  for 
phase  only  perturbations 


where  0n  is  a  real  parameter. 


4 . 1  Amplitude-only  Nulling 

Although  expensive  to  implement,  null  placement  is  possible  by  very 
few  perturbations  to  real  weights.  Adopting  a  minimax  criterion,  we 
require  to  find  the  parameters  xn  (n=1,...,N)  which  minimise  ||e§®, 
where 


Je |“  =  max  |e(u) f 
-1  <  u  <  1 

=  max  [R2(u)  +  I2(u)  )y* 

-1  s  u  £  1 

subject  to  the  constraints 

-e  S  (Re  | p ( u ) |  ,  | m | p ( u )  |  }  S  e 

for  u  =  U;  i  =  1  ,  • . . ,  p 

(where  u,  are  discrete  locations  at  which  "near-null"  placements  are 
to  be  made ) . 

The  expression  (4.4)  is  nonlinear  in  the  unknown  parameters  xn,  but 
may  be  linearized,  with  a  relative  loss  in  accuracy  of  at  most  / 2  by 
being  replaced  by 

|  e  U*  =  max  |^max|  ]  R  (u)  |  ,  |l(u)j}J 

-1  <  u  <  1 

[See  Earrodale,  Delves  and  Mason  (1978)] 

The  problem  is  now  an  overdetermined  linear  programming  problem  and 
can  be  solved  by  a  standard  routine  such  as  that  of  Roberts  and 
Barrodale  (1980). 

Figures  13  to  15  show  3  model  41  element  quiescent  patterns,  based  on 
Sine,  20  dB  Taylor  weighted  and  30  dB  Taylor  weighted  patterns, 
respectively. 

In  Figures  16  and  17  we  have  given  examples  of  null  placement  for  the 
41  element  sine  pattern.  Note  that  the  weights  are  symmetric,  thus 
reducing  the  dimensions  of  the  problem  by  half.  Also,  because  the 
perturbations  are  amplitude  only,  the  pattern  is  symmetric  about 
boresight.  In  Figure  1 6 ,  to  achieve  a  null  interval  [.7,  .73],  only 


3  pairs  o:  weights  are  changed  from  their  quiescent  values 
Figure  17  to  achieve  a  null  interval  1.07,  .08]  just  2  pairs  of 
weights  are  changed  from  unity. 

The  technique  described  here  (for  ampli tude-onlv  changes)  has  already 
been  discussed  with  different  numerical  examples,  by  Mason,  Wilde  and 
Opfer  (1987). 


9.2  phase-only  Perturbations 


ror  phase  only  perturbations,  the  constrained  pattern  is 


v  (id  u  +  e  ; 

a„  e  ■'  n 


(4.5) 


r> /  \  "  id  u  r-  id„u+  0„ 

E  (u)  =  a„  e  n  -  V  a  e  n  n 
n  n 


.  V  an  (e‘V  .  e!anu*9„, 


=  R (u)  *  i I (u) 


where 


R  ( u )  =  j  a 


n  [dos( dnu)  -  cos(dnu  *  9n) 


(4.6) 


I(u)  = 


l  an  ^sin(dnu)  -  sin(dnu  +  9n)J 


(4.7) 


Here  we  minimise 


ie  II*  =  max 


max 


(|R(u)|,  1 1  (u) | ) 


<  n .  s 


ub  s  u,  s 


for  a  discrete  set  of  u.-  (i=1,...m)  covering  the  range  of 
minimisation,  subject  to  the  constraints 


5  R(u)  S  e 


I(u)  <  e 


■>b[ua.uh]  \ 


for  a  discrete  set  of  u.-  i=m,m+1m, . . .  ,m+p  covering  the  range  o 
nulls. 

The  minimization  is  achieved  ty  imposing  inequalities 


-z  £  R  (u )  S  i. 
-z  <  Kuj)  <  z 


and  minimising  z  (in  place  of  ||e||*). 

The  inequalities  (^.9)  give 

N  "I 

-z  <  l  an  (cos(dn  ui)  -  cos(dn  ui  +  9n)}  £  z  , 

n=  1 


l 


and 


■z  <  l  an  (sin(dn  u1)  -  sin(dn  ui  +  0n))  <  z 


n=  1 


which  become: 
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:i 

:  an  cos(dn  u<)  S  I  an  cos(dn  ^  -  3n)  ♦  z  ; 

n= 1  n=1 

M  M 

l  an  cos(dn  u^)  i  _  an  cos(dn  u,  +  9n)  -  2 
n=1  n=1  ' 

(4.11) 

N  M 

l  an  sin(dn  ui)  <  l  an  sin(dn  ui;  -  5„)  *  z 
n=1  n=1  i 

N  N  ' 

I  an  din(dn  ut)  >  T  an  sin(dn  u,  +  9n)  -  2  ] 

n=1  n=  1 

Hence  the  problem  has  2p+4m  nonlinear  constraints,  ,\'  +  1  variables 
( 8n(n=1 , . . .N ) , z)  and  a  linear  objective  function  z.  ror  its  solution 
we  have  used  the  NAG  sequential  quadratic  programming  routine  E04VDF , 
as  discussed  in  §2  above. 

We  have  tackled  a  wide  variety  of  problems,  some  of  which  are  shown 
in  Figures  18-22.  The  tables  8-11  give  numerical  information  which 
corresponds,  respectively ,  to  these  5  figures. 

For  a  symmetrically  weighted  array,  if  we  limit  the  phase  changes  to 
[-ir.Tf]  then  this  results  in  conjugate  symmetric  perturbations  as  in 
the  least  squares  approach.  However,  it  is  interesting  to  note  that 
on  restricting  the  phase  perturbation  to  [0,2n],  we  produce  a  sub- 
optimal  solution  but,  by  the  nature  of  the  optimization  procedure,  few 
of  the  elements  undergo  changes.  This  is  illustrated  clearly  in 
Figures  21-22,  for  nulls  .22,  .24,  .26,  .28;  here  conjugate  symmetric 
results  are  obtained  in  Fig  21  but  all  weights  are  changed,  while 
non-symmetr ic  results  are  obtained  in  Fig  22  with  just  10  changes  to 
the  4i  weights. 


4 . 3  Loss  of  elements 


The  loss  of  elements  in  the  array  can  readily  be  counteracted  by 
applying  phase  changes  to  the  remaining  weights  so  as  to  approximate 
the  original  perturbed  pattern.  In  the  algorithms  of  §§4.1,  4.2,  we 
simply  set  the  failed  elements  to  zero  and  minimize  with  respect  to 
the  remaining  elements. 

This  technique  is  illustrated  successfully  in  Figures  23-26.  In 
Figure  23(a)  is  shown  a  perturbed  pattern  with  no  failed  element, 
obtained  from  a  20  dB  Taylor  weighted  pattern,  to  achieve  nulls  in 
[0.7,  0.72].  In  Figures  23(b)  -  23(d),  one,  two  and  three  elements 
have  failed  and  the  resulting  pattern  is  successfully  aajusted  by  the 
minimax  algorithm. 


4.4  Two-Dimensional  Arrays 

The  techniques  of  §§4.2,  4.3  are  equally  applicable  to  two- 
dimensional  arrays.  However,  the  minimization  problems  can  become 
rather  large  in  that  case,  and  so  more  efficient  but  algebraically 
complicated  techniques  such  as  that  of  Streit  (1985)  should  probably 
be  adopted  for  processing  the  linear  inequalities. 
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5. 


SUMMARY  OF  PROGRESS 


The  proposed  program  of  work  for  the  2-year  contract  was  as  follows: 

(i)  To  test  a  constrained  least  squares  method  for  adapting  a  planar 
phased  array  to  known  interference  directions 

(ii)  To  test  a  constrained  minimax  method  analogous  to  (i) 

(iii)  To  develop  and  test  new  algorithms  for  improving  (i),  (ii),  based  on 
novel  research  ideas. 

(iv)  To  extend  the  method  of  Thompson  (1976)  and  other  related  methods  for 
designing  adaptive  planar  arrays  and  to  develop  sound  nonlinear 
optimization  techniques  for  the  necessary  minimization  procedures. 

Substantial  progress  was  made  in  the  first  year  of  the  contract  on  tasks 
(i),  (ii)  and  (iii)  as  follows: 

A  modified  version,  based  on  Lagrange  multipliers,  of  Shore's  beam  space 
representation  method  for  least  squares  phased  array  adaptation  was 
introduced.  In  this  implementation  the  Lagrange  multipliers  were  in  fact 
the  beam  space  coeff icients.  The  method  was  extended  to  planar  arrays, 
and  also  successfully  applied  to  octagonal  arrays.  In  addition  multiple 
nulls  were  shown  to  be  readily  introduced,  and  formulae  for  corresponding 
Lagrange  multipliers  obtained.  Symmetry  was  considered,  and  problems  of 
pattern  enhancement  at  positions  symmetrical  to  nulls  were  studied;  it  was 
shown  that  it  is  sometimes  not  possible  to  place  symmetrical  nulls  in  beam 
space-type  algorithms. 

We  also  reported  on  work  carried  out  without  AFOSR  financial  support  on 
minimax  methods  for  null  placement  in  antenna  patterns.  We  noted  that 
algorithms,  analogous  to  those  of  Mason,  Wilde  and  Opfer  (1987)  for 
amplitude-only  nulling,  could  be  applied  to  the  phase-only  nulling 
problem.  If  phase  changes  were  restricted  to  ranges  [-ir,u]  then  the  odd 
symmetry  of  the  problem  was  adhered  to.  If  phase  changes  were  restricted 
to  [0,2ir],  then  symmetry  was  not  achieved,  but  very  few  phase  changes  were 
required  in  this  case.  It  was  also  noted  that  failed  elements  were 
readily  catered  for  in  this  type  of  approach. 

The  work  carried  out  so  far  still  requires  further  theoretical  development 
and  numerical  testing,  before  we  can  guarantee  the  full  efficacy  of  the 
techniques  discussed.  However,  the  numerical  results  produced  have  been 
consistently  good,  and  so  we  see  considerable  oromise  in  the  ideas 
introduced . 
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7  APPENDICES 


7. 1  Appendix  I 

If  it  is  assumed  that  the  phase  perturbation  will  be  small,  so  that  we  can 
use  the  linearized  form  of  the  weights  given  in  equation  (2.11),  then,  for 
symmetrical  nulls  at  ±uk ,  the  constraints  given  by  equation  (2.2)  are: 


and 


N 

j  anei<*>neldr.uk  =  0 
n  =  l 


l  ane^ne-ldnuk  =  0. 
n=1 


Putting  e^n  =  1  +  i4>n,  equations  (1.1)  become, 


•i  l  an^neldnUk  =  Po(,Jk)* 

n=1 


(1.1) 


and  (1.2) 

N 

~i  l  an<one-ldnuk  =  p0(uk)# 
n=»1 

Here  pQ(uk)  is  the  value  of  the  unperturbed  pattern  at  uk.  Clearly  the 
left  hand  sides  of  equations  (1.2)  are  complex  conjugates  and  due  to  the 
odd  symmetry  of  4>n  and  d^  the  resulting  equations  are  inconsistent. 
Therefore  it  is  not  possible  to  gain  a  solution  using  the  linearized  form 
for  symmetric  null  placement. 


7.2  Appendix  II 


Given , 


F  =  2  I  cn|anl 


:os> 


n=  i 


then  the  derivative  is  given  by, 


—  =  2  l  cn|an|2sinon  lin  ,  p=1 , 
3Ap  n-1  3Ap 


(II. 1 ) 


Also , 


3C.  „  V  3'i>n  s  =  1  ,  .  .  P , 

3Ap  n=  1  3 <p n  3Ap  P-1...P. 


(II. 2) 


and  itzs  is  given  by  equations  (2.21)  and  (2.22). 

3dn 

From  6 QUow  ion  (2.23)  v»g  **ovg  9 


^n  = 


tan 


-I  1  Xkv(dn)Vsin[dn(un'us)]  *  l  2  Xkv{dn)Vcos[dn(’Jk 
k=l  v=0 , 2 . .  k*1  v=1 , 3- • 

step  2  step 

2cn I an I  *  l  l  Akv(dn)Vcos[dn(uk-us)]  "  l  l  Xkv(dn)Vs 
k=1  v=0,2..  k=l  v“1,3«. 

(II. 3) 


Putting 

then 


Yn  =  tan 


lin  =  9y. 


aXp  dyn  3AP 


(II. 4) 


(II. 5) 


and 


■^n  =  cos2i(i 
3yn 


(II. 6) 


Denoting, 


vn  "  2=nlanl 


+  l  l  Xkv(dn)V«s[dn(uk-us)]  -l  l  Akv  (dn)^sin[dn( 

k  even  v  k  odd  v 


and 


tn  -  -l  l  Xkv(dn)vsinCdn(uk-u3)]  *  I  l  Xkv(dn)'Jcos[dn(uk-us)l, 
k  even  v  k  odd  v 


^  “  -  (in)Vsin[dn(uk-us)]  +  (dn)vcos[dn(ak-us)]  In 


for  even  v,  (II. 7) 


and 

(ln)Vcos[dn(uk-us)]  -  (dn)vsin[dn(uk-us)]  in 

V  v 

n  n  for  odd  v  (II. 8) 

Equations  (II. 7)  and  (II. 8)  together  with  (II. 6)  (2.21),  (2.22),  (II. 2) 
and  (II. 1)  give  ail  the  required  derivatives  for  the  objective  function 
and  constraints. 


oooo  nano  on 


7.3  Appendix  III  Listing  of  Routine  for  Setting  non-Octagon  weight 
to  Zero. 


SUBROUTINE  WEIGHT { A ,  C,  NELEM,  NADD,  NSOCT ,  MAXELT ,  POLY) 

C 

C  IN  THIS  ROUTINE  THE  WEIGHTS  A  AND  C  ARE  SET  ACCORDING  TO 

C  THE  TYPE  OF  ARRAY,  RECTANGULAR  (R)  OR  OCTAGONAL  (O). 

C  NELEM  IS  TOTAL  NUMBER  OF  ELEMENTS  IN  RECTANGULAR  ARRAY. 

C  NSOCT  IS  NUMBER  OF  ELEMENTS  ALONG  EDGE  OF  OCTAGONAL  ARRAY. 

C  NADD  IS  NUMBER  OF  ELEMENTS  ADDED  TO  EACH  SIDE  OF  OCTAGONAL 

EDGE  TO  EMBED  IT  INTO  RECTANGLE. 

IMPLICIT  DOUBLE  PRECISION  ( A-H ,  P-Z) 

CHARACTER* 1  POLY 

DIMENSION  A(MAXELT) ,C (MAXELT) 

FIRST  SET  ALL  WEIGHTS  TO  1,  FOR  SINC  PATTERN 

DO  10  JELEM  =  1,  NELEM 
C(JELEM)  =  l.OdO 
A (JELEM)  =  l.OdO 
10  CONTINUE 

IF  ( LIT . EQ . ’ R * . OR . LIT . EQ . ’ r ’ )  RETURN 

NOW  PLACE  ZEROS  FOR  OCTAGONAL 

DO  20  JELEM  =  1 ,  NADD 
A (JELEM)  =  0.0D0 
A(NELEM-JELEM+1 )  =  0 . 0D0 
20  CONTINUE 
C 

NSTART  =  NSOCT  +NADD 
NUMAD  =  NSOCT 
NZERO  =  2*NADD  -1 
C 

DO  25  INUM  =  1,  NADD 
DO  23  JELEM  =  1,  NZERO 
NTJ  =  NSTART  +  JELEM 
A(NTJ)  =  0.0D0 
A( NELEM-NTJ+ 1 )  =  0 . 0D0 
23  CONTINUE 

C 

NUMAD  =  NUMAD+2 
NSTART  =  NTJ+NUMAD 
NZERO  =  NZERO  -2 
25  CONTINUE 
C 

RETURN 

END 


List  of  Figures 

Figure 

1  Geometry  of  1-D  array 

2  Enhancement  of  pattern  of  Symmetric  Location  to  null: 
u=0. 43633 

3  Placement  of  symmetric  nulls  at  locations  between  nulls 
of  quiescent  pattern 

4  Increasing  null  width  by  higher  order  constraints 

5  Coordinate  system  for  2-D  array 

6  Elements  of  2-D  array 

7  Embedded  octagon  in  2-D  array 

8(a)  Field  pattern  for  rectangular  array:  13x13  elements 
8(b)  Field  pattern  for  octagonal  array  embedded  in  13x13 
rectangular  array 

9(a)  Null  placed  at  u=0.28,  v=0.32 

9(b)  Perturbed  -  Quiescent  pattern  for  rectangular  array 
10(a)  Null  placed  at  u=0.28,  v=0.32 
10(b)  Perturbed  -  Quiescent  pattern  for  octagonal  array 
11(a)  Nulls  placed  at  u=0.28,  v=0.32  and  u=0.32,  v=0.36 
11(b)  Perturbed  -  Quiescent  pattern  for  rectangular  array 
12(a)  Nulls  placed  at  u=0.28,  v=0.32  and  u=0.32.  v=0.36 
12(b)  Perturbed  -  Quiescent  pattern  for  octagonal  array 

13  Quiescent  sine  pattern  with  41  elements 

14  Quiescent  20  db  Taylor  pattern  with  41  elements  and  n=6 

15  Quiescent  30  db  Taylor  pattern  with  41  elements  and  n=6 

16  Constrained  amplitude-only  pattern.  Nulls  at  u=0.7,0.71. 
0.72,0.73 

17  Constrained  amplitude-only  pattern.  Nulls  at  u=0.07,0.08 

18  Constrained  phase-only  pattern.  Nulls  at  u=0.07,0.08 
Phase  range  [-11,11] 

19  Constrained  phase-only  pattern.  Nulls  at  u=0.4,0.525, 
0.65,0.775,0.9.  Phase  range  [-11,11] 

Constrained  phase-only  pattern.  Nulls  at  u=0.22,0.24. 
0.26,0.28.0.3,0.32,0.34.0.36.  Phase  range  [-11.111 


20 


21  Constrained  phase-only  pattern.  Nulls  at  u=0.22,0.24. 
0.26,0.28.  Phase  range  f-IT ,IIJ 

22  Constrained  phase-only  pattern.  Nulls  at  u  =  0.22.0.24, 
0.26,0.28.  Phase  range  (0,211] 

23(a)  Constrained  phase-only  pattern.  Constraints  over 
u=[0.7,0.72] 

23(b)  Array  element  4  failed 
23(c)  Array  elements  4,11  failed 
23(d)  Array  elements  4,11,21  failed 


Fig.  5 


Coordinate  system  for  2-D  array 


Fig.  6 


Fig.  7 


Elements  of  2-D  array 


Embedded  octagon  in  2-D  array 


■  -  -  - . 

» j.  !  ~  2 1**0  -  -2v 


•  f  ♦  •  f-  T  ■ 


2(u 


T 


Fig.  12(a) 


Nulls  placed  at  u=0.28.v-0.32  and 


u=0.32.v=0.36 


Fig.  12(b) 


Perturbed -Quiescent  pattern  for  Octagonal  Array 


Jm. 


Constrained  amplitude -only  pattern. 

Nulls  at  u=0.7,0. 1 1 ,0. 1 2,0.  <  3 

Quiescent  pattern  -  sine 

Xi=  1.000000000  (i=l,...,41)  except 

Xi=i4i=0.27835309 

x2— 2^40— 0.65933325 

x3  =  £39 =0.96323350 


Fig.  16 


XI  0  1 


Fig.  IT  Constrained  amplitude-only  pattern. 
Nulls  at  u=0.07,0.08 
Quiescent  pattern  -  sine 
Xi  =  1.000000000  (i=l,...,41)  except 
x5=ai37=-0. 61903885 


List  of  Tables 


Table 

1(a)  Results  for  null  at  u=0. 43633 

2(a)  Results  obtained  for  multiple  nulls 

2(b)  Effect  of  close  null  locations 

3  Coefficients  and  average  perturbations  for  higher 
order  constraints 

4  Results  for  rectangular  array  with  one  null 

5  Results  for  octagonal  array  with  one  null 

6  Results  for  rectangular  array  with  two  nulls 

7  Results  for  octagonal  array  with  two  nulls 

8  Effect  of  2  near  boresight  constraints  on  30  db 
Taylor  pattern.  Phase  range  [-11.11]  (Fig.  18) 

9  Effect  of  5  widely  spaced  constraints  on  30  db 
Taylor  pattern.  Phase  range  [-11.111  (Fig.  19) 

10  Effect  of  8  closely  spaced  constraints  on  30  db 
Taylor  pattern.  Phase  range  [-11. II]  (Fig.  20) 

li(a)  Effect  of  4  constraints  on  30  db  Taylor  pattern 
Phase  range  [-11.11]  (Fig.  21) 

11(b)  Effect  of  4  constraints  on  30  db  Taylor  pattern 
Phase  range  [0,211]  (Fig.  22) 


Mu 1 1  direction  Optimal  Beam  Overage  Phase  Null  Depth  Loss  at  Main 

(rad)  UoBl'f icients  Pertubation  (rac)  (DBR)  Beam 


Null  direction  Optimal  Beam  Average  Phase  Null  Depth  Loss  at  Main 

(rad)  Coefficients  Pertubation  (rad).  (DBA)  Beam 


Table  4-  Results  for  Rectangular  Array  with  one  Null. 


Antenna  Arrays  :  Phase  Onl£  Nulling 

There  are  169  elements  and  1  contraints. 

Constraint  1  is  at  position  .28000  .32000 

Initial  value  for  constraint  1  is  .10000 

Required  accuracy  is  . 10000E-07 

Optimal  Beam  Coefficients 
. 15315E-01 


Coefficients  for  Perturbed  pattern 


1 

. 72655E-02 

2 

. 59292E-02 

3 

- . 967  10E-03 

4 

- . 695  14E-02 

5 

- . 6439  IE-02 

6 

- . 48876E-17 

7 

. 64391E-02 

8 

. 69514E-02 

9 

. 967 10E-03 

10 

- . 59292E-02 

1  1 

- . 72655E-02 

12 

- . 18904E-02 

13 

. 52129E-02 

14 

. 6492  IE-02 

15 

. 57311E-17 

16 

- . 6492  IE-02 

61  . 18904E-02 

17 

- . 69062E-02 

62  . 72655E-02 

18 

- . 95252E-03 

63  .  59292E-02 

19 

. 58716E-02 

64  -.967 10E-03 

20 

.73000E-02 

65  -.6951 4E-02 

21 

. 191 86E-02 

66  -  .  44733E-02 

22 

- . 527 14E-02 

67  .  27 990E-02 

23 

-.75111E-02 

68  .7511  IE-02 

24 

-.27990E-02 

69  .5271 4E-02 

25 

. 44733E-02 

70  -.19 1 86E-02 

26 

.7646  IE-02 

71  - . 7  3000E-02 

27 

. 967 10E-03 

72  -  .  587  16E-02 

28 

- . 59292E-02 

73  . 95252E-03 

29 

-.72655E-02 

74  . 69062E-02 

30 

-. 18904E-02 

75  .6492  IE-02 

31 

. 521 29E-02 

76  -  .  7 6805E-18 

32 

•75327E-02 

77  -.649  2  IE-02 

33 

.28392E-02 

78  -  .  69062E-02 

34 

-.45291E-02 

79  . 18904E-02 

35 

-.76387E-02 

80  . 72655E-02 

36 

- . 36645E-02 

81  . 59292E-02 

37 

.  36645E-02 

82  -.9671 0E-03 

38 

.  7 6 387E-02 

83  - . 695  14E-02 

39 

. 45291E-02 

84  -.64391E-02 

40 

- . 52714E-02 

85  .00000 

41 

-.7511  IE-02 

42 

- . 27990E-02 

Average  pertubation  =  .48752E-02 

43 

.447  3  3E-02 

44 

.  7646  IE-02 

45 

. 37 140E-02 

Quisecent  pattern  at  uk  =  .28000 

46 

- . 37 140E-02 

Real  =  -.64964  Imag  —  . 

47 

- . 76461E-02 

48 

- . 44733E-02 

Perturbed  pattern  at  uk  =  .28000 

49 

. 27990E-02 

Real  =  .  1 2688E-09  Imag  = 

50 

. 75111E-02 

51 

. 527 14E-02 

Null  depth  =  -240.14  DBA 

52 

- . 19186E-02 

53 

- . 76387E-02 

54 

- . 36645E-02 

55 

.  366  45E-02 

56 

. 76387E-02 

57 

. 4529  IE-02 

58 

- . 28392E-02 

59 

- . 75327E-02 

60 

- . 52 129E-02 

51033E- 
1186 1E- 


1 


1 


.  32000 
.32000 


Table  S  Results  for  Octagonal  Array  with  one  Null. 
Antenna  Arrays  : 

There  are  129  elements  and  1  contraints. 

Constraint  I  is  at  position  -28000 

Initial  value  for  constraint  1  is  .10000 

Required  accuracy  is  .10000E-07 

Optimal  Beam  Coefficients 

.27348 


32000 


Coefficients  for  Perturbed  pattern 
1  O  0  Q  A 


5  .12394 

6  .10187E-15 

7  -.12394 

8  -.11639 

q  -  .  I  5089E-0 1 

17  .13062 

18  .  1  9825E-0 1 

19  -.11491 

20  -.12413 

21  - . 300 1 9E-0 1 

22  .8491 4E-0 1 

23  .13698 

29  .13496 

30  .3917  6E-0 1 

31  - . 10359 

32  -.13022 

33  _  .  4 4629E-0 1 

34  .72241E-01 

35  .13679 

36  .7  47  02E-0 1 

37  -  .  747 02E-0  1 

41  .13698 

42  . 57604E-01 

43  -  .  90124E-0  1 

44  -.13449 

45  -  .  58757E-01 

46  . 58757E-01 

47  .13449 

48  .901 2  4E-0 1 

49  - . 57604E-0  1 

50  -.13698 

51  -.84914E-01 

53  .13679 

54  .  7 4702E-0 1 

55  -  .  74 70  2E-0 1 

56  -.13679 

57  -  .  7224 1E-0 1 

58  . 44629E-0  1 

59  .13022 

60  .10359 

61  -.391 76E-0 1 

62  -.13496 

63  -.9661 OE-O 1 

64  .  1  5089E-0 1 

65  .11639 

66  .901 24E-0  1 

67  - . 57604E-01 

68  -.13698 

69  - . 8  4  9 1 4E-0  1 

70  .3001  9E-0  1 

71  .12413 

72  .11491 

73  -  .  1 9825E-0 1 

74  -.13062 

75  -  .  10716 

76  .  11972E-16 

77  .10716 

78  .13062 

79  - . 39 1 76E-0 1 

80  -.13496 


81 

82 

83 

84 

85 


.  966  10E-0 1 
.  15089E-01 
.11639 
.12394 
.00000 


Average  pertubation  =  .67559E-01 

Quisecent  Pa^^6at  uk  imag  =  '  28°°°2 547 9E-  1 6 

Perturbed  patterngat =  ^ 

Null  depth  =  -342.29 


.28000 

.  1415  IE-15 


DBA 


.  32000 
. 32000 


Table  6 


Results  for  Rectangular  Array  with  2  Nulls 


Antenna  Arravs  :  Phase  Only  Nulling 

There  are  169  elements  and  contraints  ■ 


Constraint 

Constraint 


1  is  at  position 

2  is  at  position 


.  28000 
.32000 


Initial  value  for  constraint 
Initial  value  for  constraint 


1  is 

2  is 


. 10000 
.  10000 


Required  accuracy  is 
Optimal  Beam  Coefficients 
.  86483E-0  1 
-  .  89406E-0 1 


. 1 0000E-0' 


Coefficients  for 


1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 

14 
1  5 
16 

17 

18 

19 

20 
21 
22 

23 

24 

25 

26 

27 

28 

29 

30 

31 

32 

33 

34 

35 

36 

37 

38 

39 

40 

41 

42 

43 

44 

45 

46 

47 

48 

49 

50 

51 

52 

53 

54 

55 

56 

57 


53802E- 
1 1934E- 
.  46955E- 
38105E- 
37794E- 
,  3393  IE 
,  25908E- 
.  49  172E- 
.  20995E 
.  1 1646E- 
.  34537E- 
.  57008E 
.  10019E- 
.62571E- 
.  43776E 
.  41216E- 
.  13255E 
.  31713E 
.  28 18  IE 
.  14293E- 
.20031E- 
.  13013E- 
.21166E- 
.56119E- 
. 86034E- 
.75122E- 
.  39847E- 
.43719E- 
.64120E- 
. 29024E- 
. 29987E- 
. 20832E- 
. 18744E- 
. 14180E- 
.74719E- 
. 54348E- 
.70521E- 
. 14511E- 
.89014E- 
. 45568E- 
.  11409E- 
. 25907E- 
.  3  1  302E- 
.  55600E- 
■ .  17150E 
■ . 15 130E 
•. 6 3 36  IE 
. 51725E 
■ .  53892E 
■ . 21284E 
. 80033E 
.  156  10E 
• .  16247E 
.  22414E 
.  32112E 
.  894  18E 
15275E 


Perturbed  pattern 
01 


02 
0  1 
01 
02 
0  1 
01 
02 
01 
01 
02 
02 
02 
02 
01 
01 
02 
01 
01 
02 
01 
01 
02 
02 
03 
03 
01 
01 
02 
01 
01 
02 
01 
01 
03 
02 
03 
02 
02 
01 
01 
01 
-01 
-02 
-01 
-01 
-03 
02 
03 
02 
02 
-01 
01 
-01 
-01 
-02 
-01 


58  15846E-01 

59  -  .  20047E-02 

60  .  4829  IE-02 

61  -  .  36409E-03 

62  -  .  27727E-02 

S3  .69819E-02 

64  .16  17  5E-0 1 

65  .  5608  IE-02 

66  .  18598E-0 1 

67  . 32407E-0 1 

68  . 1 2  1 7 1E-0  1 

69  -. 13147E-01 

70  -.1631 5E-0  1 

71  - . 33453E-02 

72  . 44 1 03E-02 

73  - . 18350E-03 

74  - .  337  39E-02 

75  . 58532E-02 

76  .  16483E-0I 

77  .82735E-02 

78  -  .  17961E-01 

79  .  32191E-01 

80  .  1  5194E-01 

81  -  .  10800E-0 1 

82  - .  16529E-01 

83  - . 46347E-02 

84  . 39226E-02 

85  .00000 

Average  pertubation  = 


. 32000 
.  36000 


,  14419E-01 


Quisecent  pattern  at  uk 
Real  =  -.64964 

Perturbed  pattern  at  uk  = 

Real  =  -.191 10E-07  Imag  = 


-  . 28000 

Imag  =  .  51033E-15 


.  28000 

.  15856E- 


15 


Null  depth  = 


-196 . 59 


DBA 


Quisecent  pattern  at  uk  -  .32000 

Real  =  .84424  Imag  = 


34094E-15 


Perturbed  pattern  at  uk  - 
Real  =  - .  20194E-07  Imag  = 


.  32000 

- . 12334E-14 


Null  depth  = 


-196.11 


DBA 


.  32000 
.  32000 

. 36000 
. 36000 


Table  7  Results  for  Octagonal  Array  with  two  Nulls. 


Antenna  Arrays  :  Phase  Onl  £  Nulling 

There  are  129  elements  and  2  contraints. 

Constraint  1  is  at  position  .28000  .32000 

Constraint  2  is  at  position  .32000  .36000 

Initial  value  for  constraint  1  is  .10000 

Initial  value  for  constraint  2  is  .10000 

Required  accuracy  is  .  10000E-07 

Optimal  Beam  Coefficients 
- . 59158 
. 37090 


Coefficients  for  Perturbed  pattern 

5  .10705 

6  -.17203 

7  -.19683 

8  - . 7  3508E-0  1 


9 

. 62827E-0 1 

17 

. 14239 

18 

- . 14722 

19 

- . 20278 

20 

-.9110  1E-0  1 

21 

. 46339E-0 1 

22 

.12795 

23 

. 11257 

29 

. 17337 

30 

-. 11836 

31 

- . 20582 

32 

- . 10779 

33 

. 29  386E-0 1 

34 

.11997 

35 

.11885 

36 

.  37  580E-0  1 

37 

- . 58779E- 0  1 

41 

. 19944 

42 

- . 86054E-0  1 

43 

-.20569 

44 

-. 12337 

45 

. 12141E-01 

46 

. 11046 

47 

.12300 

48 

. 52316E-01 

49 

- . 45232E-01 

50 

-.  10985 

51 

-.10702 

53 

. 22033 

54 

-.51 102E-0 1 

5  5 

-.20213 

56 

-  .13764 

57 

- . 52235E-02 

58 

. 99605E-0  1 

59 

.  12505 

60 

. 65926E-0  1 

61 

- . 30714E-01 

62 

-. 10472 

63 

-.11400 

64 

- . 46184E-01 

65 

. 72443E-01 

66 

- . 14530E-01 

67 

-. 19490 

68 

- . 15037 

69 

- . 22537E-01 

70 

. 87  57  8E-0  1 

71 

.  12508 

72 

. 78154E-01 

73 

- . 15528E-01 

74 

- . 97679E-0 1 

75 

-. 11943 

76 

- . 60697E-01 

77 

. 56321E-0  1 

78 

. 17034 

79  -.18386 

80  -.16135 

81  - . 39627E-0  1 

82  . 7  4  5  5 1E-0 1 

83  .12316 

84  . 88790E-01 

85  .00000 

Average  pertubation  = 


Quisecent  pattern  at  uk 
Real  =  9.0426 

Perturbed  pattern  at  uk 
Real  =  -  .  46  4  1  3E- 1 1 

Null  depth  =  -268.8 


. 80078E-0  1 

=  .28000 

Imag  =  - . 25479E-16 

=  .28000 

Imag  =  -.46632E-15 

DBA 


Quisecent  pattern  at  uk  =  .32000 

Real  =  4.7164  Imag  =  -.11891E-14 

Perturbed  pattern  at  uk  =  .32000 

Real  =  -.58462E-12  Imag  =  .46453E-15 

Null  depth  =  -286.87  DBA 


. 3200C 
. 3200C 

. 3600C 
. 3600( 


Table  8 


Effect  of  2  near  boresight  constraints  on  30  db 
Taylor  pattern.  Phase  range  [-11,11]  (Fig.  18) 


BEAM  REFERENCE  AT  CENTRE  OF  ARRAT 
INITIAL  PATTERN  TATL3R  WEIGHTED 
NUMB  cR  UP  EC  UAL  SlOcLOBES 
DEPTH  of  S I0EL36ES 
NUMBER  of  elements 


NUMBER  QF  FITTING  POINTS 
OEPTM  OF  NULL 
NULL  INTERVAL 
SCALING  FACTOR 
E“SI(NULL  OEPTM  *  SCALING 

CT0LIA6SJLUT  E) 
FTOLCRELATIVE) 

CONSTRAINT  NUMBER 

1 

2 

INTERVAL  OF  Phase  CmangES 


6 

30.0 
4  I 
4  1 

-90. 

CO. 070,0. 0803 
I  .  000 
FACTOR);  *  592TE-33 
=  . 1000E-05 
;  .  100DE-06 

VALUE  TMETA 

0.07 00  -.01 

0.0800  4.59 


objective  function 

norm  Cp  TmE  constraint  VIOLATIONS 


2  .  5  7  7  A  A  I  3  ♦  3  0 
1.2902190-08 


NO. 


OP  ITERATIONS  =  19 

cpu  time  csecs) 

OBJECTIVE  FUNCTION 
GA  IN 

number  of  phase  perturiat IONS 
AVERAGE  PMASc  PERTURBATION 


1  0  •> .  9  3  0 
2.677 

-0.3399 

h  1 

0.00000030 


[  WEIGHTS 

1  0.264375373671 

2  0.272106558834 

J  0.258425718517 

4  0.3  1448  10  52  305 

5  0.350356583175 

6  0.3  96  3979  3  1  568 

7  0.450536777323 

3  0.508724360485 

9  0 . 563 25-6 15 102 

10  0.626581310598 

H  0 .632261419286 

12  0.734354271079 

13  0.73441 2668475 

1 4  0.830863725639 

11  0.3TJ62 166 8153 

16  0 .9  1  16009o9  0  4  8 

17  0.943565035459 

13  0.968566193776 

19  0.936207431334 

20  0.996588198007 

21  1 . OOOOCOOOOOOO 

2  2  0.9  96533  1  98007 

2  3  0 . 9  3  620  74 J 1 88  4 

2  4  0.9  66  56  61  98  776 

25  0.943565035459 

26  0.911600969043 

2  7  0 .8  7362  16681  58 

28  0.830863725839 

29  0.734412668475 

30  0. 73435 4Z 71079 

3  1  0 .682  26  14  19  286 

32  0.626581310593 

33  0.5o8254615102 

34  0.503724360485 

35  0.450536777328 

3  6  0.3  9689  79  31  666 

37  0.350356533176 

38  0.314481052303 

39  0.238425713517 

4  o  0.2  72  1065  59334 

41  0.264375373671 


PHASE  CHANGES 
-0.458410754363 
-0.249570927309 
-0.373585951 613E-01 
0.164306088971 
0.35074  2503  194 

0.513156764071 

0.63616 1908346 
0.693381875330 
0 .640  51  07  5  1  069 
0.44038  3273  382 
0. 13011 3553453 
0.114J65339513E-01 
-0  .4  3599  92  24  330E-0  L 
-0 .3499  1  6940  58  3  6-0 l 
-0. 288275961 065E-02 
0.29565972774 9 £-0 l 
0.501396197S81E-01 
0 .5  4404  1965  62  6E-0 1 
0.44251 5323909E-C 1 
0.243283195082E-01 
-0.1  7817  7  3  3454  36-1  5 
-0.2432831 95082E-0 1 
-0.44251 5323809E-01 
-0.54494 1955626E-0 1 
-0 . 50139  61  97  38  1  E-0 1 
-0 .29565  9727  7496-0 l 
0.238175961065E-02 
0. 349916940533E-01 
0.435 999224 830E-01 
-0.114066339518E-G 1 
-0.130113553453 
-0.440383273382 
-0.640510751 069 
-0.69338 1875330 
-0 .63616 1908846 
-0.513156764071 
-0. 350  74  2  50 J  1  94 
-0.164306088971 
0.3  7  353  59  5  1  6  1  3  E-0 1 
0.249570927307 
0.468410754369 


F 


7* 


Table  9 


Effect  of  5  widely  spaced  constraints  on  30  db 
Taylor  pattern.  Phase  range  [-11,11]  (Fig.  19) 


3£AN  REFERENCE  47  CENTRE  Q  F  ARRAT 

Initial  pattern  taylor  weighted 

NUNS  £3  OP  EQUAL  SI0EL0SE3  =  6 

OEPTH  UP  SIOELDBES  --  30.; 

number  op  elements  a  4i 

NUMBER  OP  PITTING  POINTS  -  41 

OEPTH  OP  NULL  =  -00. 

NULL  INTERVAL  =  CO. 400, 

SCAL ING  FACTOR  =  1 .000 

epsicnull  depth  *  scaling  factors 

CTOLCAaSQLUTE) 
ptolcrelativeo 
Constraint  number 
i 


0.0000 


Interval  op  phase  changes 


*  c 


Exit  eoavcp  -  optimal  solution  pound. 

objective  FUNCTION 
norm  of  the  constraint  violations 

NO.  QP  ITERATIONS  ■=  3  3 

CPU  time  OECSl 
OBJECTIVE  FUNCTION 
GAIN 

NUMBER  CF  PHASE  PERTURBATIONS 
average  phase  perturbation 


5927  =-G3 
10005—05 
10005-06 

V4LU5 

THETA 

0.4000 

23.53 

0 . 5250 

3  1.67 

0.5500 

40.54 

0-7750 

50.3  1 

0.9000 

-pi  >  ?n 

6-.  16 

1.2132170 

-01 

1  .  20  56490 

-0  3 

*  525-900 

=  .1212 

=  -0.9398 

=  4  V 

0 . 00000000 


I  WEIGHTS 

1  0. CS4375373&T1 

2  0.27210*559534 

3  0.288425710517 

4  Q. 314  lb  1052  30  5 

5  0.35085*533176 

6  0.396397931*63 

7  0 .45Q536TT7 323 

8  Q. 5087243*0485 

9  0 . 5o825»6 15 10 2 

10  0.626581810593 

11  0-682261419286 

12  0.734J54271079 

13  0-784412663475 

14  0.830863725339 

15  0.873621663153 

16  0 .9 11 nO 096904 J 

17  0.943565035450 

13  0.9*856619377 5 

19  0.936207431384 

20  0.9965831  53007 

21  1.000000000000 

22  0.996588138007 

23  0.936207431334 

24  0.96856*1  18  77* 

25  0.94356303*459 

28  0.911600939043 

2  7  0.8  7362  16  88  153 

23  0.830363725839 

29  0 . 7444  12863475 

30  0.734354271079 

31  0.632261419288 

32  0.626531810518 

33  0 . 5682 546 15 102 

34  0. 5087243*0485 

35  0.45053*777323 

36  0.396997931863 

37  0.350358533176 

38  0.314481052305 

39  0.233426718517 

40  0.272106559834 

41  0.264375373671 


PHASE  CHANGES 
0.692330413065E-D1 
0.20209  5009  21  1 
-0  .  87511  2684  18  7  E-0  1 
-0 .4*607  0249  28  5E-0 1 
-0  .  1*697  14  2707  3 
0 .20948  99648  39 
-0.333589401 507E-01 
-0.900932690993E-01 
0 .3  3602  1793  505E-0 l 
-0 .677396290  47  3E-0 l 
0 .12239  1362  484 
-0 . 858*4 1794006E-0 I 
0.362373398464E-01 
-0-705272236369E-01 
-0.3771360360T7E-02 
0.143176434634 
-0.137735778942 
0-253295522052E-0 1 
-0.41  13  8  36  37492E-0 1 
0.102125135364 
0.32810 13092046-19 
-0.102  12  5  1  3  5  24  ) 

0  .411)3  J629452E-0 1 
-0-253295522446E-01 
0.  13T7  3  5  T78119 
-0.14317*404339 
0.3771  3  59  1043  5E-02 
0  .  TO  5  2  7  22  34  2O5E-01 
-0.3623783986276-01 
0.838641794508E-01 

-0.1  2239  1364  74Q 

0.6773962938065-0 1 
-0.33602 1794186E-01 
0 .900  13  26897  32E-0 1 
0.383389401 307E-01 
-0 .209483984640 
0 .16697  1427  09  1 
0 .46607  024941  7E-0  1 
0.8  751  1  2634  553E-01 
-0.202395007367 
-0.6323304090636-01 


Table  10  Effect  of  8  closely  spaced  constraints  on  30  db 
Taylor  pattern.  Phase  range  [-11,11]  (Fig.  20) 


beam  REFERENCE  at  centre  c  F  array 
initial  pattern  taylcs  -eighteo 


NUMd  6  R  OF  E2U&L  SI0SL38ES 
0 EPTH  OF  5I3cL03£i 

NUMBER  OF  ELS ME  MTS 

NUMBER  OF  FITTING  POINTS 
QrPTH  OF  NULL 

NULL  INTERVAL 

-  6 
=  30.0 

=  Hi 

-  HI 

-  -90  . 

»  CO. 220.0. 

360  : 

scaling  factor 
fc  PS  I  ( NULL  DEPTH  *  SCALING 
C  fOLC ABSOLUTS ) 
FTQLCRELATive) 

constraint  .NUMBER 

«  1.000 

F1CTCR)=  .5927E-03 
-  .  10006-05 
=  .10006-02 
value 

THETA 

1 

0.2200 

12.71 

2 

0 . 2  h  0  0 

13.83 

2 

0.2600 

15.07 

* 

0.2900 

16.26 

s 

0.3000 

17.46 

6 

0.3200 

13.66 

7 

0 . 3  HOO 

10-33 

8 

0 . 3  600 

2  1.10 

interval  of  phase  changes  =  c-pi,pi3 

exit  cqv/of  -  current  point  cannot  be  improve 

3  UPON. 

objective  function 

NORM  OF  TrtE  CONSTRAINT  VIOLATIONS 

NO.  CF  ITERATIONS  =  6- 

C=U  TIME  (SECS) 

OBJECTIVE  FUNCTION 
GAIN 

NUMJER  OF  PHASE  PERTURBATIONS 
AVERAGE  PHASE  perturbation 


a  .  9  35  59-D-O 1 
1 .532*  3  70-  1  5 


397.710 
.  393  5 

-0  .2151 
*  1 

0 . 0  0  0  9  0  0  0  C 


I  WEIGHTS 

1  0.26*375373671 

2  0.272106553834 

3  0.233*26713517 

*  0.31-431052305 

5  0.350856583176 

6  0.396897931668 

7  0.450536777328 

8  0.50872*360*35 

9  0.56825*615102 

10  0.626531310598 

11  0.632261-19286 

12  0.73*85*271079 

13  0.78**12663*75 

1*  0.830363725839 

I5‘  0.37  362  16681  58 

16  0.9116009690*3 

17  0.9*3565035*59 

13  0.963566193776 

19  0.986207*3188* 

20  0.996538198007 

21  1.000000000000 

22  0. 996538198007 

23  0.986207-3138* 

2*  0.963566193776 

25  0.9*3565035*59 

26  0.9116009690*8 

27  0.873621663158 

28  0.830863725839 

29  0.78**12668*75 

30  0. 73*854271079 

31  0.6822oi*19286 

32  0.626531310593 

3  3  0  .  5682  5*6  1  5  10  2 

3*  0.50372*360*35 

35  0.450536777328 

36  0.396897931668 

37  0.350856583176 

38  0.31**31052305 

39  0.238*26719517 

*0  0.27210655983* 

*1  0.26-375373671 


phase  Changes 

1  .  71  5  10  27  7802 
0  .  1  719  56  5  1  500  7  E-01 
0  .  3662  02  *9  *8  36 
0.2*1*  33*79*  66 
-0.1*32  3*3**925 
-0.103902811326 
0  .  18*3  15  2667  62E-0 1 
0  .  36  9*  2*67  7*  27  5-  0  1 
0  .  13*6*368  3360E-0 1 
-0. 611 1 3 1* 1 225 3E-0 1 
-0  .  11  17  5*98  7377E-0 1 
0.  *3 06  53318237  5-  01 
0.*90l 56101 261E-0 1 
-0.46-1 35535-125-02 
-0  .*39*  6562  93  985-0  1 
0  .  160  1  18  7*  38165-0  2 
0  .  3296  53  55  3  1  325-  0  1 
0.37  20  23  21  1  1  7  SE-02 
-0.213953676*385-01 
-0.239*05317080E-01 
-U  .  1  526  J137  625  8E- 1 0 
0.239*053181155-01 
0.2139536758675-01 
-  0  .  3720  23  19  72  28  E-02 
-0  .  3296  53  35  5  7325-  0  1 
-0  .  160  1  1  8  7*  3  7  3*  5-0  2 
0.439*656295125-01 
0.46-1 35536530E-02 
-0.490156101*005-01 
-0.4356583137965-01 
0.11176*986*155-01 
0.61113 1* 12*675-0 1 
-0 . 1 3*6-86832695-01 
-0.389*2*6774005-01 
-0.  18*3  1  5  2667  7  35-0  1 
0.103902311332 
0.1*823*3**823 
-0.2*1*33*79515 
-0. 365202*9*8*3 
-0.171956515138E-01 
-1.71510278011 


r1  . . .  * 

f 


►  Table  11(a)  Effect  of  4  constraints  on  30  db  Tavior  pattern 

[  Phase  range  [-11,11]  (Fig.  21) 


t 

► 


L 


SEAM  REFERENCE  AT  CENTRE  QF  IRJiY 
INITIAL  PATTERN  TAYLOR  WEIGHTED 


NUMBER  OF  equal  s; 

C6L0BE  S 

6 

OEPTH  OF  SIDcLQBES 

= 

30  •  0 

number  of  elements 

- 

4| 

NUMBER  QF  FITTING 

POINTS 

- 

41 

NUMBER  OF  C3NSTRAI 

NTS 

- 

4 

OEPTH  QF  NULl 

- 

-90. 

NULL  INTERVAL 

- 

CO. 220,0* 

230  3 

SCALING  FACTOR 

- 

1  .000 

EPSICNULL  DEPTH  x 

SCALING 

FACTOR  )  = 

. 5927E-03 

CTDLCABSOLUTE) 

= 

.  1000  E-05 

FTOLCRELAT IVE) 

= 

.  1  000  =-06 

CONSTRAINT  NUMBER 

VALUE 

TnE  7  A 

I 

0 .2200 

12.71 

2 

0.2400 

13.39 

3 

0.2600 

15.97 

*» 

0. 2300 

16.26 

INTERVAL  Or  PHASE 

CHANGE  S 

C -P  C , PI  o 

08J6CT:v£  FUNCTION  '  5 . 3029 6 5C-0 1 

NORM  OF  Trt£  CONSTRAINT  VIOLATIONS  =  3.^*^  5  5  3  30-  16 


NO.  OF  ITERATIONS  =  330 

CPU  TIME  C  SECS  3 
OBJECTIVE  FUNCTION 
GAIN 

NUM8  c  R  OF  PHASE  ? £ R T UR3 A T I  0 N S 
AVERAGE  ?*ASE  PERT1JR3ATI  ON 


1  5  7  2  .  720 
.  5803 

-0.1534 
4  1 

0. 00030000 


\ 


► 

r 


l 


t 


i  u;:c.hts 

1  0.25-375370671 

2  0 .2 72  1  0  65  59  83- 

3  0.293-26718517 

-  0.31--31052305 

5  0.350355533176 

6  0.396397931663 

7  0.-50536777323 

3  0 . 508  72-350-d  5 

9  0.56825-615102 

10  0.626581510593 

11  0 .632261-19236 

12  0.73-35-271079 

13  0.78--12658-75 

1-  0.830363725339 

15  0.873621663153 

16  0.9116009690-3 

17  0.9-3565035-59 

13  0.963566138776 

19  0.956207-3133- 

20  0.996583193007 

21  1.000000000000 

22  0.996533199007 

23  0.936207-3138- 

2-  0.9o8S66196776 

25  0.9-3565035-59 

26  0.9116009690-3 

2  7  0 . 3  7  36  2  l  663  3  53 

23  0.330863725339 

2  3  0  .73--1  2663-7  5 

30  0.73-35-271079 

31  0.632261-19286 

32  0.626581310599 

33  0.56825-615102 

3-  0.50872-350-35 

35  0.450536777323 

36  0.396397931663 

37  0.350353533176 

33  0.314-81052305 

39  0.238-26718517 

-0  0.27210655933- 

-1  0.26-373373671 


PM-S5  C-UN053 
1.-236 1892--3 
0.1  7730  533--8  3 
0.5---110-3385E-01 
0 .3523803507195-0  1 
-0.2171333389535-01 
0.  16008  53  754835-0  2 
-0 . 3  7731  3-  5-  508  5-0  1 
-0  .  3  323-  27  20  22  15-0  1 
0 . 20768-3600055-0 1 
0.  13023  075-  3705-02 
-0.708933298-735-02 
-0.3-45650832635-03 
-0-102599228359E-01 
0. 2  0  8  6 1  9030  -356-0  2 
0.36-7151707576-02 
-0.-333669152215-02 
0 .9-652  6631  7--E-0  3 
0 .2  5002  -106  19  25-0  2 
-0.-286336-0 1975-02 
0.2573856022195-02 
0 .29-3  7  1  3  5  1  7725-  1  - 
-0.2575356022205-02 
0  .-  2363  36-0  19  7E-02 
-0.2  5002  -1  06  19  25-0  2 
-0 .9-632  6691  7-06-0  3 
0  . -3386  59  1  52215-02 
-0.3647151707536-02 
-0.208619030-365-02 
0  .  102699226  359E-0 1 
0.3*4  56  5083  26  56-0  3 
0 . 70393  3298-755-02 
-0.18028075-3705-02 
-0.2076843600055-01 
0  .  3  32  3-  2  720  22  15-0  1 
0.37731 3-5-6035-0 1 
-0.160085375-796-02 
0.2171333889525-01 
-0.3523803507206-01 
-0  .  5-4-1  10-3  3875-0  l 
-0.17730583--86 
-  1  .-256  18  92--3 


ft 


Table  11(b)  Effect  of  4  constraints  on  30  db  Taylor  pattern 
Phase  range  [0,211]  (Fig.  22) 
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i 

i 

) 

[ 
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Q F  4  R  R  A  r 
WEIGHTED 


BEAM  REFERENCE  a  f  CENTRE 
INITIAL  PATTERN  T-TLJR 
NUMBER  OF  cwUAL  SIJELQ3ES 
DEPTH  OF  SIDEL08lS 

number  of  elements 

NUMBER  OF  FITTING  POINTS 
NUMBER  of  constraints 
depth  of  null 

NULL  INTERVAL 

scaling  factor 

e?5I(NULl  06PTH  X  SCALING 
CTOLCABSOLUT  z) 
FT3LCRELATIVS) 

CONSTRAINT  number 
1 
z 

3 

interval  of  phase  Changes 


30.0 
A  i 
: 


2  . 

=  c  3.2: 3, 0- 2S0 J 
-  1 .000 


F AC  TOR  )  =  .  5  92  7  5-03 
=  .10005-05 
=  . 1000E-06 
VAL  J€ 

THETA 

0.2200 

12.71 

0.2400 

13.99 

0.2600 

15.07 

0. 2B00 

1  o.  26 

*  CO.OOPI3 

OBJECTIVE  FUNCTION  =  l .  3- 37Q60* 30 

NORM  OF  THE  CONSTRAINT  VIOLATIONS  =  1.419:320-03 


NO.  CF  ITERATIONS  =  23 

CPU  TIME  (SECS) 

QoJECTIVE  FUNCTION 
GA  IN 

NUMBER  CF  PHASE  PERT'JRSA  TI3NS 
average  PHASE  PERTURBATION 


73.930 

1  . 

*0.4411 

i  a 

0 . U  $  0  1  0  7  2 


► 


r 


I 

WEIGHTS 

1 

0.2  69  3  7  5  3  7  36  7  1 

o 

0.272106559939 

3 

0.283920 718517 

*♦ 

0.319491052305 

5 

0.350356533174 

6 

0.396397931669 

7 

0.950536777323 

8 

0  .  508  729360985 

9 

0 . 5  632594  1  5  1  02 

10 

0.626531310593 

1  1 

0 .632  26  19  19  286 

12 

0.739359271079 

13 

0.739912668975 

1  4 

0.830363725339 

1  5 

0.873421668153 

1  6 

0.911400969098 

1  7 

0.993565035959 

13 

0-963566198776 

19 

0.936207931389 

20 

0.996593198007 

2  1 

1 .000000000000 

22 

0.996583198007 

23 

0.936207931389 

24 

0.968566193776 

23 

0.943565035959 

26 

0.911600969093 

27 

0 .8  7362  1668  1  53 

2  3 

0.830363725339 

29 

0 .739.12648975 

30 

0.739359271079 

3  1 

0.632241919286 

32 

0 .62653 1810598 

33 

0.548259615102 

34 

0.508729340935 

35 

0.950536777328 

36 

0 .396897931668 

3  r 

0.350356533176 

38 

0.314931032305 

39 

0.238426718517 

40 

0.272106559334 

4  1 

0.264375373671 

P  M  A  S ;  CHANGES 
j  .  1  6  3  1  i  O  5  5  3  J 
0.647  50  95  3  3  0  72 
1  . 563  2  5  9*5 454 
0.t2C)l 1443415 
0.0000000000006*00 
0.0000000000006*00 
0.0000000000006*00 
0.0000000000006*00 
0.0000000000006*00 
0.0000000000006*00 
0.0000000000006*00 
0.0000000000006*00 
0.0000000000006*00 
0.191025312276 
0.9285709975216-01 
0.0000000000006*00 
0.0  0000  0000  0006  *00 
0.0000000000006*00 
0.0000000000006*00 
0.0000000000006*00 
0.0000000000006*00 
0.0000000000006*00 
0.0000000000006*00 
0.0000000000006*00 
0.0000000000006*00 
0  .  1  4  “•  3  1  0  5  3  6  3  5  9 
0.569  59  93  96  79  56-0  1 
0.0000000000006*00 
0.0000000000006*00 
0.  00000  00000006  *00 
0. 00000  00000006  *00 
0.0000000000006*00 
0.0000000000006*00 
0. 183ol 8278172 
0.0000000000006*00 
0.0000000000006*00 
0.0000000000006*00 
0.0000000000006*00 
2.01159099199 
0 . 00000  0000  0006*00 
0.0000000000006*00 
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